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Abstract 

We investigate in a relativistic Hartree theory the gross properties of cold 
symmetric nuclear matter and nuclear matter in beta equilibrium under the 
influence of strong magnetic fields. If the field strengths are above the critical 
values for electrons and protons, the respective phase spaces are strongly 
modified. This results in additional binding of the systems with distinctively 
softer equations of state compared to the field free cases. For magnetic field 
~ 10 20 Gauss and beyond, the nuclear matter in beta equilibrium practically 
converts into a stable proton rich matter. 
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In recent years considerable efforts have been directed to the study of the effects of intense 
magnetic fields on various astrophysical phenomena. Large magnetic fields B m = 10 12 — 10 14 
G have been associated with the surfaces of supernovae |1[ and neutron stars On the 

other hand, extremely large fields could exit in the interior of a star. It is presumed from the 
scalar virial theorem || that the interior field in neutron stars could be as high as ~ 10 18 
G. Besides, the matter density in the neutron star core could exceed up to a few times the 
nuclear matter density. At such high fields and/or matter density, constituents of matter are 
relativistic. Moreover, the energy of a charged particle changes significantly in the quantum 
limit if the magnetic field is comparable to or above a critical value. The critical field is 
defined as that value where the cyclotron quantum is equal to or above the rest energy 
of the charged particle, which for electrons is B^ c > = 4.4 x 10 13 G, and for protons it is 
Bffi^ ~ 10 20 G. Theoretical studies of free electron gas in intense magnetic fields relevant to 
the neutron star crust have been carried out by several authors using the Dirac theory || as 
well as Thomas- Fermi and Thomas- Fermi-Dirac models [H. The intense fields were shown 
to drastically reduce photon opacities and greatly accelerate the cooling rates in neutron 
stars ||. It has been also demonstrated |7],|] that the magnetic fields have significant effects 
on the weak interaction rates and the abundances of light elements in the early Universe. 
The influence of extremely large fields on neutron matter || relevant to the neutron star 
interior and on the thermodynamic properties of strange quark matter in cosmic QCD phase 
transition and baryon inhomogeneity in the early Universe |10j have been also investigated. 

Motivated by the existence of strong magnetic fields which quantize the motion of the 
electrons, we investigate in this Letter its influence on the gross properties of dense nuclear 
matter appropriate to the interior of a neutron star. This may have profound implications 
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on cooling rates, mass-radius relationship of neutron stars. It is also instructive to extend 
the calculations to values of B m > 10 20 G where along with the electron, the proton motion 
is strongly quantized. Fields of such magnitude, appropriate to neutron star interior, could 
largely modify the proton phase space in the quantum limit. Though such high field is 
hitherto unestimated, it may possibly exist in the core of neutron star. 

We therefore consider strong magnetic field effects on nuclear matter and a system com- 
posed of neutrons, protons and electrons (n-p-e system) in beta equilibrium within a rela- 



tivistic Hartree approach in the linear a-co-p model [11]. In the beta equilibrium case, the 



electrons are assumed to move freely in the strong magnetic field, whereas the produced 
neutrinos/anti- neutrinos escape from the system without being Pauli blocked. In a uniform 
magnetic field B m along z-axis, the relativistic Hartree Lagrangian is given by 



C = tp 



hn D ^ ~m- g a a- - ^g P l^r ■ p M 



k=w,p 



M " l ml{Vif 



(1) 



in the usual notation [11]. Here, _D M = <9 M + iqA^, where the choice of gauge corresponding 
to the constant B m along z-axis is A Q = 0, A = (0, xB m , 0). The general solution for protons 
is ip(r) oc e~ leHt+tPvV+vpzZ f Py , Pz {x) , where f PytPz {x) is the 4-component spinor solution. The 
Dirac-Hartree equation for protons in a magnetic field is then given by 



-ia x d/dx + a y (p y - qB m x) + a z p z + (3m* + U" p \ f^ y ] Pz (x) = e H f^ Pz (x 



'Py,Pz 



(2) 



The equation of motion for neutrons is obtained by setting the charge q = and replacing 
U^ p by U^ n in Eqs. (1) and (2); the corresponding solution is a plane wave. It may be 
mentioned that the Dirac theory for free electrons in a homogeneous magnetic field was first 
studied by Rabi fll2| , and can be obtained by putting U$ = in Eq. (2). Since we confine 
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to cold systems (T = 0), only positive energy spinors are considered. These in the chiral 
representation are of the forms 



( (e? +p z )I u . Py (x) \ 
-iy/2vqB m I u - 1]Py (x] 



-m*I v . Py (x) 



(3) 



-m*I„- lxPy (x 



-iy/2uqB m I v . Py {x] 
\(e? +p z )I u _ 1;Py (x)J 



(4) 



where N v = l/^/2e^(e^ +p g ), and = e H - U$ p = (p 2 z + m* 2 + 2uqB m ) 1 / 2 is the effective 
Hartree energy, with v the Landau principal quantum number which can take all possible 
positive integer values including zero. The function I u - Py (x) is similar in form as in Ref. 
|T3| . The effective nucleon mass m* = m + U$, where the nucleon rest mass is taken 



m p = 939 MeV, and Ug = —{g a /m a ) 2 n s . The total scalar density is 



as m = m n 



n s = n s + n s'\ with 



n 
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2^ 
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m*qB m ^ x 

2_, 9u m 



2tt 2 



u=0 



U* + 01/2 



{m* 2 + 2vqB m ) 1/2 



(5) 



(6) 



where O n = //* — m* , and O p , v = fi* — m* — 2vqB m . The interaction energy density 
Uq for protons and neutrons are given by U^ p = (g u /Tn u ) 2 TiB + {.9p/ m p) 2 P^/^ an d U^ n = 



{gu/m w ) 2 n B - (g p /m p ) 2 p 3 /A, where p 3 

n B — n n + With 



n n . The total baryon number density is 



n r . 



03/2 
n n 



r< JJ "max 

— - V a 1/2 

o„2 y^ P ,u ■ 

Z7r i/=0 



(7) 



Here z/ max is the largest integer not exceeding (/x* 2 — m* 2 )/(2qB m ), and the effective chemical 
potential /z* is at the Fermi surface. The Landau level degeneracy factor g v is 1 for v = 
and 2 for v > 0. The total energy density of the system is given by 
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Here = p? e — m 2 — 2uqB m and m* 2 = m* 2 + 2quB m , where m*s denote m*s(m e s) for 
i = p(e). The first three terms of Eq. (8) correspond to the interaction energy densities 
for er, and p mesons. The last three terms are the expressions for kinetic energy densities 
for neutrons, protons and electrons. The total pressure generated by the system is given 
by P = n 2 B d(E/A)/dn B , where E/A is the total energy per baryon. For symmetric nuclear 
matter (where n n = n p = n^/2), m* is evaluated self-consistently for a given % and 
B m . On the other hand, the n-p-e system under the beta equilibrium and the charge 
neutrality conditions is in particular important for neutron star. For these two cases, when 
B m > B$(°\ the charge neutrality condition, n p = n e , gives 



. (e) 



E 9»oy 2 = e 9,oi; 2 . 



(9) 



u=0 



u=0 



When B m > B$(°\ but appreciably smaller than B^°\ large number of Landau levels are 
populated and the relations are almost similar to the field-free case. However, when B m 
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significantly affects the electrons so that z/ max is small (« 0), the protons are also affected. 
(An estimate of z/^ x f° r various values of B m is discussed later in the text.) Employing 
Eq. (9) in conjunction with the /3-equilibrium condition, \i n = \i v + /i e , one can obtain m* 
self-consistently for a given ub and B m . The proton and neutron chemical potentials, fi p and 
pL n , are related to their respective fermi momenta k^ and k^ by pL p = +[kp^ + m l 2 ^] 1 ^ 2 
and //„ = Ufr n + [k^ + m* 2 ] 1//2 . Therefore, the neutral p meson field affects the chemical 
composition inside the neutron star through the different proton and neutron vector potential 
Uq in the asymmetric n-p-e system. 

In the present calculation the parameters for the coupling constants and mesons masses 
are taken from Horowitz and Serot ]14| to be ^(m/m ff ) 2 = 357.47, g 2 (m/m u ) 2 = 273.87, 
and g 2 (m/m p ) 2 = 97.00. This yields nuclear matter saturation density at no = 0.1484 
fm~ 3 with a binding energy of 15.75 MeV and a bulk symmetry energy of 35 MeV. In the 
top panel of Fig. 1, the variation of effective nucleon mass m*/m with baryon density 
ub/uq is displayed. The curves (a) and (b) represent the symmetric nuclear matter case 
for B m = and B^", respectively. It is found that for B m = 0, m* decreases gradually 
with ub while for B^ c ', the decrease is relatively much faster beyond % m n . This is 
attributed to the drastic reduction in the proton fermi momentum k^ , whereas the neutron 
fermi momentum k^ is unaffected by B m and is identical to the fermi momentum, kp for 
B m = 0. Consequently, at any tib, p* p is smaller than /i*. These are reflected in the larger 
value of ns and hence in the magnitude of Ug for non-zero magnetic field, in contrast to the 
field free case. By further increasing B m to 1QB^ C \ m* jm for symmetric nuclear matter 
(curve (c)) undergoes a further reduction beyond the density ~ 3no. 

Considering now a n-p-e system, the curve (d) in the top panel of Fig. 1 shows the 



variation of m* jm for such a system at B m = 0. If B m < B^°\ the variation of m* 
remains virtually unaltered from the field free case (not shown in the figure). If the field 
is further increased to B^^ and 10B^°\ the m*/m values are significantly reduced as 
evident from curves (e) and (f) of the figure; the distinction between them occurs only at 
ub > 2n . Furthermore, in presence of B m , the m* values here are found to be much smaller 
than those for symmetric nuclear matter. This may be attributed to the neutron-proton 
asymmetry in the system. 

The remarkable distinction of m* and kp^ in intense magnetic field from those of the 
field free case is expected to be manifested in the equation of state (EOS) which is crucial 
in understanding the gross properties of dense matter. The energy per baryon E/A with 
varying baryon density n B /n is exhibited in Fig. 2. The curves (a) and (b) respectively, 
correspond to B m = and lOB&l^ for symmetric nuclear matter. It is observed that the 
strong magnetic field B rn > B$^ causes additional binding of the nuclear matter with 
E/A ~ -41 MeV for B rn = 10B^ C \ The usual binding energy curve for n-p-e system in 
absence of magnetic field (curve (c)) shows no binding. When the field is slightly quantizing, 
the system is still unbound as indicated by curve (d) for B m = 10 4 B^; it is only sightly 
softer than (c). However, when both protons and electrons are strongly quantized by B m , 
the n-p-e system is strongly bound, and the binding increases with B m as exhibited by 
curves (e) and (f) for B rn = B$^ and 10B^ C \ respectively. In contrast to B m = 
case, for non-zero B m , even though the contribution from the scalar density is increased, the 
relatively larger decrease in kinetic energy density especially for protons results in the excess 
binding. Furthermore, it is observed that with increasing B m , the minima of the binding 
energy curves, where the pressure P = 0, shift towards higher densities. This is clearly seen 
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in the insert of Fig. 2 where the pressure P is displayed as a function of energy density e; the 
curves (a) to (f) correspond to the same values of B m as in Fig. 2. The causality condition 
dP/de < 1 is fulfilled by all the cases considered here. It is evident from Eq. (8) that 
the kinetic energy density for protons is strongly suppressed, and a meson term is strongly 
enhanced in the magnetic field. The latter term has a negative contribution to the pressure. 
On the other hand, u and p meson terms (p3 = for symmetric matter) which increase 
with tib, compensate the reduction in the kinetic energy and the scalar meson terms in the 
pressure at higher density to produce zero pressure (or energy minimum) compared to the 
B m = cases. For the n-p-e system, considerable suppression of kp and m* in magnetic 
field accentuates the above effect, as a consequence it is more bound with the minimum 
occurring at a higher density than the symmetric nuclear matter. 

Recent studies have indicated that proton fraction in neutron star matter is crucial in 
determining the direct URCA process which leads to the cooling of neutron stars |T5|,[TE|]. In 



the bottom panel of Fig. 1, the proton fraction Y p = n p /riB is shown for the n-p-e system 
for B rn = (solid line) and for lO^Bffiw (dashed line). The proton fraction is observed to 
be enhanced in the latter case. For direct URCA process, the inequality kp + kp > kp \ 



which corresponds to Y p > 0.11 for B m = HlBj , should be satisfied. In the linear a-uo-p 



model with B m = 0, this condition is satisfied at n B > 1.5n and thus rapid cooling by 
direct URCA process can occur. On the other hand, for B m = BwW, the proton fraction 
shown by the dotted line in the figure, is found to be considerably enhanced. The drastic 
fall in the proton fermi momentum entails a substantial n — > p conversion, as a result the 
system is converted to a highly proton rich matter. Moreover, it has been demonstrated in 
Fig. 2 that such systems are energetically more favorable. Therefore, if the magnetic field 
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is strong enough ~ 10 20 G, possible existence of stable "proton matter" may be envisaged. 
If the field is further increased to 10B^ C \ the proton fraction (shown by dash-dotted line) 
saturates to a value of 0.98 at ub > 2 no- 

When B m > B^ c \ v^l x for various values of ub and B m is found to follow the relation- 
ship v^l x ~ [1 / \B m j B^^)}[X(nB /no) — J7(n£/n ) 2 ], where for symmetric nuclear matter, 
X = 101217.12 and J = 5458.64, and for the n-p-e system in beta equilibrium X = 46571.24 
and J = 562.35. Thus for a fixed decreases monotonically with increasing B m . For 

all ub values of interest, = when _B m ^10 6 B^ c \ and this is found to be in conformity 
with the values B m > B^^ (/i e /m e ) 2 /2 predicted in Ref. Ul7| . As a consequence of charge 



neutrality, is^l x = v^lx = an d ^f = ^f ( see Eq. (9) and Ref. Therefore in such 
strong fields, the direct URCA process in stars would occur if Y v > (X 1 / 3 — X^ 1 ^ 3 /3d), which 
corresponds to the real positive root of the above mentioned inequality condition. In this 
expression X = {l + (l + 4/27rf) 1 / 2 }/2rf, and d = U^n B /2>{qB m f . Interestingly, Y p depends 
not only onng but also on B m . Comparing the values of the proton fraction obtained from 
the model calculations (Fig. 1) and the inequality condition, we find the threshold for direct 
URCA process is not reached for B m > B^ c \ 

The effect of intense fields on the neutron star profiles is obtained by applying the 
EOS to solve the Tolman-Oppenheimer-Volkoff equation [|ll|. For magnetic fields B m = 0, 



lO^WM, B^ c \ and 105 <W( C >, the maximum masses of the stars are found to be M max = 
3.1OM , 2.99M Q , 2.91M , and 2.86M Q , respectively. The corresponding radii are -RM max = 
15.02, 14.95, 12.25, and 12.00 km. These values suggest that the neutron stars masses are 
practically insensitive to the effects of the magnetic fields, whereas the radii decrease in 
intense fields, leading to their compactness. 
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In this letter, we primarily focus on the new qualitative features that arise out of nuclear 
matter in a strong magnetic field within a relativistic Hartree approach in a simple linear o- 
uj-p model. We believe that these features will survive even in more sophisticated calculations 
with a more refined EOS. It will be worth investigating the influence of a quantizing field 
on the quark matter in a relativistic Hartree-Fock model. 
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Figure Captions 



FIG. 1. The variations of effective nucleoli mass m*/m (top panel) and proton fraction Y p 
(bottom panel) with baryon density n B /n for different values of magnetic field B m discussed 
in the text. 

FIG. 2. The energy per baryon E/A as a function of n B /n for different values of B m . 
In the insert, the pressure P is shown as a function of energy density e for different values 
of B m (for details, see text). 
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